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Kibble mechanism drastically underestimates the production of topological defects, as confirmed recently in 
atomic and condensed matter systems. If non-thermally produced, they can be cosmological dark matter of mass 
1-10 PeV. If thermalized, skyrmion of mass 1-10 TeV is also a viable dark matter candidate, whose decay may 
explain e~ spectra in cosmic rays recently measured by PAMELA, FERMI, and HESS. Models that produce 
magnetic monopoles below the inflation scale, such as Pati-Salam unification, are excluded. 


Topological defects are of common interest to condensed 
matter physics, atomic physics, astrophysics and cosmology, 
as well as algebraic topology (il). When the symmetry group 
G spontaneously breaks down to its subgroup H, there are 
continuously connected ground states parametrized by the 
coset space G/H. The homotopy groups of the coset space 
then tell us what kinds of topological effects are possible. 
In most cases, non-trivial 7q(G/H) implies the existence of 
(2 — d)-dimensional topological defect. If the coset space 
has disconnected pieces (7(G/H) #4 0), we expect domain 
walls. For multiply-connected space (7(G/H) 4 0), there 
are strings (vortices). If the boundary of space can map non- 
trivially to the coset space (72(G/H) # 0), we expect point- 
like defects such as magnetic monopoles. An exception to the 
rule is when the whole space is mapped non-trivially to the 
coset space (773(G/H) # 0), where skyrmions are stabilized 
by non-renormalizable terms in the low-energy effective the- 
ory (3. In this case, it is not the boundary condition that is 
topologically non-trivial, but the configuration in the bulk. 


To estimate the initial abundance of defects produced by a 
phase transition in early universe, Kibble pointed out that the 
correlation length diverges at the critical temperature while 
the causality does not permit exchange of information beyond 
the horizon scale He therefore came up with a lower 
bound on the amount of defects, namely approximately one 
per horizon, called Kibble mechanism. Most of the litera- 
ture uses this lower bound as the estimate of the abundance 
of topological defects from phase transitions in early uni- 
verse. For point-like topological defects, one finds nrp/s ~ 
(T../Mp1)?. Therefore, only phase transitions close to the 
grand-unification scale produce abundance of topological de- 
fects worthy of consideration. 


A decade later, Zurek [5|] proposed a more refined estimate 
of the abundance by carefully considering the time scale avail- 
able. His estimate has been confirmed experimentally in a 
large number of systems recently, now called Kibble—Zurek 
mechanism. The studies include liquid crystals , super- 
fluid *He igi and ‘He ALG, an optical Kerr medium [11], 
Josephson junctions 12411 3h), superconducting films (14, and 
spinor Bose-Einstein condensate 15). 


We point out that the Kibble—Zurek mechanism provides 
a substantially larger abundance of topological defects from 


phase transitions in early universe than the original estimate 
by Kibble. Therefore even phase transitions just above the 
TeV energy scale may produce interesting (or dangerous) 
amount of topological defects. 

In particular, we discuss the possibility that point-like topo- 
logical defects may be the cosmological dark matter, which 
is arguably one of the most pressing mysteries in cosmol- 
ogy, astrophysics, and particle physics i BA. The domi- 
nant paradigm to explain the dark matter is the thermal relic 
of yet-undiscovered particle. Within this paradigm, we con- 
sider dark matter candidates below approximately 100 TeV in 
mass because of the unitarity bound . Our main result in 
this Letter is that the natural range for topological dark matter, 
if non-thermally produced by a second-order phase transition, 
is O(1 ~ 10) PeV, which obviously violate the unitarity limit. 
Note that a symmetry breaking at this energy scale in the hid- 
den sector is of great interest in many attempts to understand 
the origin of hierarchy between the Planck and electroweak 
scales such as dynamical supersymmetry breaking, and extra 
dimensions. In addition, we also point out that skyrmions at 
the order 10 TeV, once thermalized, are also interesting dark 
matter candidates that are often ignored in the literature (171. 
The existence of skyrmion solution is quite generic in mod- 
els where Higgs serves as a pseudo Nambu Goldstone boson, 
which opens the new possibility to connect the origin of elec- 
troweak symmetry breaking and dark matter. 

If the dark matter particles are produced thermally at tem- 
peratures higher than their mass, their initial abundance is the 
same as any other relativistic particle species. Then the final 
abundance is determined by their annihilation cross section, 

1.1 x 10°(¢ + ie, Cay 3x 10~?7cm3/sec 
Q X hk? y e AS Be a a X — 
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a) 
where xf = m / Ty with Tẹ the freeze-out temperature, and 
we used g, œ~ 100 and £ = 0 (S-wave). Assuming that only 
one partial wave J would contribute, the annihilation cross 
section is limited from above by 


4n(2J +1) 3 x 107? (2J + 1)cm?/sec 
(m/TeV)? 


Combining Eqs. UB), we find m < 110 TeV assuming S- 
wave annihilation and J = 0. 
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On the other hand, the Kibble—Zurek mechanism predicts 
a very different abundance of point-like topological defects. 
Throughout this paper, we assume second-order phase transi- 
tion. The correlation length € and relaxation time 7 diverge 
near the critical temperature which can be parametrized using 
the critical exponents v and u respectively 


E ~ Eole” , 


where € = (Te — T)/T, characterizes the proximity to the 
critical temperature Te. 

The system is quenched when it passes through the criti- 
cal temperature with a finite speed. It is characterized by the 
quenching rate tg = (t — t.)/e to the linear order around 
time te when T = Te. During the quenching, there exists a 
particular time t, when the time remaining before the transi- 
tion equals the equilibrium relaxation time |t, — te| = T(t). 
Beyond this point the system can no longer adjust fast enough 
to follow the changing temperature of the bath, and at time 
t,. the fluctuation becomes frozen until a time |t,. — te| af- 
ter the critical temperature is reached. It is easy to see that 
le(tx)| = (te/7)~1/" +"). Therefore, the fluctuation does 
not get smoothed out beyond the correlation length 


Elta) ~ £0(TQ/T0) TF . (4) 


1/2 


T= Tle“; (3) 


In radiation dominated universe, T «x t7 and one finds 
TQ = 2te = H(T.)~' with the expansion rate H = a/a. 
Assuming the free energy of the Landau—Ginzburg form 
V(¢) = (T — Temo? + $A¢* near Te, one can approxi- 
mate m ~ XT, and £, 7 scale as £0/\/€, 70 /,/€ classically. 
So the critical exponents are u = v = 3. Setting the initial 
correlation length ĉo © To ~ 1/(VATe), we have 


1/3 1 2 1/3 
Ex Te ——_ = H! H (5) 
H ALT, AT? 


In radiation dominated universe 


T? 45 
H = — = 
CMp ’ G V 413g. (6) 


and hence the correlation is shorter than the horizon size by 
a factor ~ (T;/Mp;)?/°, leading to a far larger number of 
defects than the original Kibble’s estimate. 

For point-like defects (PD), we expect approximately one 
per £3. Assuming (g, œ 107 — 103, A ~ 0.3-1), we find 


Te 
PED a0, 
Mpi 


(7) 
s IT=T, 

If we consider the quantum corrections to the system, the 
critical exponents u and v could be different from 1/2. Gen- 
erally speaking, it is related with the anomalous dimension of 
the leading relevant operators in the of the Lagrangian of the 
scalars that triggers the symmetry breaking. Causality € < cr 
dictates v < u and we will assume v = u below as the Hubble 
friction term for scalar field ¢ could be ignored in the vicinity 


2 


of critical temperature indicated in footnote (40). For typi- 
cal quantum systems based on O(N )-symmetric ¢* theory in 
three dimensions fil, the critical exponents are v = 0.625 
in binary liquid system (N = 1), v = 0.672 in superfluid 
tHe experiment (N = 2), and v = 0.70 in EuO, EuS system 
(N = 3). As we can see, v is quite close to 2/3 and does 
not vary very much for different N. By plugging in the same 
numbers as the classical case, we obtain 


(8) 


We see that the magnetic monopoles are produced orders of 
magnitude more than the original Kibble’s estimate and hence 
even models with phase transitions down to TeV scale are sub- 
ject to serious constraints. For instance, Pati-Salam model 
|36] assumes the symmetry breaking SU (4)o x SU(2)z x 
SU(2)r —> SU(3)o x SU(2)z x U(1)y, and hence pre- 
dicts magnetic monopoles. Once produced, the monopoles 
are stable and their number can only be reduced by annihila- 
tion of M — M pairs. The magnetic (hidden) monopoles will 
stay in kinetic equilibrium with the (hidden) thermal plasma 
of electrically charged particles. The long-range forces be- 
tween M and light charged particles will dissipate the energy 
of M drifting towards a nearby M, allowing capture and an- 
nihilation to occur. Preskill found the annihilation was 
negligible assuming the initial abundance given by the Kibble 
mechanism. With the Kibble—Zurek mechanism, however, the 
annihilation must be considered for magnetic monopoles, and 


we must use Eq. (5) in Ed), 


nM 1 4n\? mpp _992 Te 
= = 7.9 x 107? 
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where B = (3/477)¢(3) >, (hqi/4r)? which sums over all 
spin states of relativistic charged particles and h = 27/q is 
the magnetic coupling. 

In Fig. [1] combing Eq. (8) and (9), we show how magnetic 
monopole density to entropy ratio depends on the critical ex- 
ponent v in 2nd order phase transition for different critical 
temperatures. It is clearly that the Parker limit excludes 
the such monopoles assuming the phase transition below the 
unification scale unless the critical exponent v is significantly 
above 1. However, this is not the case normally considered in 
relativistic field theories for phase transitions. 

On the other hand, the point-like defects may be magnetic 
monopole under a U(1) gauge theory unrelated to electro- 
magnetism (“hidden U(1)”). We assume that there is no cor- 
responding hidden plasma to dissipate the energy of M and 
we ignore the annihilation. We can also ignore the annihi- 
lation for point-like defects based on global symmetries be- 
cause there is no long-range force among them. For non- 
thermal production of topological defects to dominate, they 
have to be heavy enough so that they never stay in chemical 
equilibrium once produced. At the critical temperature, when 
Le = Mpp/T. > xf, which ranges from 20 to 30 for differ- 
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FIG. 1: nm/s versus the critical exponent v if the magnetic 
monopoles are produced by a 2nd order phase transition. The dashed 
line are the naz /s without the annihilation. The green region is al- 
lowed by the Parker’s bound naz/s < 10~?° while the blue region 
is allowed by the direct search nm/s <2x 10778 (38). 
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FIG. 2: Relic density of topological dark matter as a function of 
its mass based on Eq. if it is non-thermally produced during a 
second order phase transition. The yellow horizontal band denotes 
the relic density 0.094 < Qh? < 0.129 preferred by WMAP data. 
We assume ze = mpp/Te = 50. 


ent monopole coupling strengths and phase transition temper- 
atures, the relics density could be derived from Eq. as 


3v 

T.\ (30T.\ #7 
Opph? = 1.5 x 10° ( Ze £ . 10 
PD i (= Mp a 


If we take xe = 50, the relic density is a function of Te, as 
shown in Fig.|2| In order to account for the cold dark matter 
abundance, we need T, ~ O(1) PeV to in the classical case 
and T, ~ O(10) PeV in the typical quantum cases (42). 
What kind of model can lead to realistic topological dark 
matter? One obvious possibility is that there is a hidden non- 


TABLE I: Summary of popular composite Higgs models in 3 + 1 
dimensions that generate skyrmions. 


Models L<« [4] 13(G/H) 
Minimal Moose [20] [su su) | Z 
Littlest Higgs [21] | su) |506) | Zo 
SO(5) Moose [22] | s06)?] 506) | Z 


abelian gauge theory whose breaking to U(1) produces mag- 
netic monopoles. As long as the U (1) does not mix with QED, 
strong bounds such as Parker’s limit do not apply. Their 
annihilation cross section in the plasma is negligible. As an 
example, the vector-like model of dynamical supersymme- 
try breaking by Izawa—Yanagida and Intriligator-Thomas 

has SO(6) global symmetry. Gauging SO(3) subgroup 
embedded diagonally into SO(3) x SO(3) c SO(6), one can 
see that it breaks to SO(2) and produces magnetic monopoles. 

It is not clear if skyrmions can be created by the same 
Kibble—Zurek mechanism, as they are topologically non- 
trivial configurations in the bulk rather than the boundary con- 
ditions. However, skyrmions are baryonic composites in the 
underlying gauge theory and hence may be thermalized inde- 
pendent of the production mechanism. Note that many com- 
posite Higgs models (e.g., little Higgs theories) proposed in 
the literature can have skyrmions as topological solitons (see 
Table (43]). Their masses are expected in the 10 TeV re- 
gion because fy ~ 1 TeV in these theories from the natural- 
ness argument. Once thermally produced, the correct abun- 
dance of topological dark matter could be obtained with a rel- 
atively strong coupling gpp ~ 3. Since the global symme- 
try G in those models in Table [lis approximate, we may ask 
whether the skyrmion is metastable. Gauging a subgroup of 
G may induce the skyrmion to decay through instanton effects 
(271. However, the enormous suppression factor proportional 
to exp(—87?/g) will make its life time much longer than the 
age of our universe (44). 

Let us now comment on the consequence of topologi- 
cal dark matter on cosmic ray signals. For the case of 
skyrmion dark matter, we can imagine that skyrmions will de- 
cay through some higher dimension operators analogous to 
proton decay in Grand Unified Theory (GUT). The most eco- 
nomical way is to consider GUT-suppressed dimension 6 op- 
erators, with its lifetime 


M4 TeV \° 
re Sr TU = 3 x 107s ( x) ( 
mpp MPD 


Maur i 
2 x 10!6GeV 
aul 


The final decay products of the skyrmions would be some 
meson states with extra fundamental fermions, for instance 
charged leptons. We can imagine that the main branching ra- 
tio of the skyrmion decay is the one into a light meson state 
below GeV. The light meson mixes with the Higgs boson, so 
its coupling to the SM fermions is proportional to their masses 
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and will dominantly decays into u pairs. As long as the mass 
of the skyrmion is multi-TeV, the muon dominated leptonic 
final state will naturally explain the PAMELA excess 
in et while no excess in p [32] as well as slow decline 
E~*-° in the e~ + e+ spectrum as reported by FERMI [33] 
which steepens at about | TeV as measured by H.E.S.S. [34]. 

In summary, we have considered the possibilities that point- 
like topological defects, such as monopoles and skyrmions, as 
the viable dark matter candidates. We apply the Kibble—Zurek 
mechanism to the non-thermal production of monopoles by 
a second order phase transition, and find that the abundance 
is much larger than the one originally estimated by Kibble. 
Depending on critical exponent in the correlation length, the 
hidden monopoles could account for the correct relics den- 
sity for the mass range of approximately 1—10 PeV. The ther- 
mally produced skyrmion of mass 1—10 TeV can also provide 
the correct relics density of cold dark matter, whose decay 
may account for the flux and spectral shape of electrons and 
positrons recent observed by PAMELA and FERMI satellites 
and HESS. A similar consideration excludes any models that 
produce magnetic monopoles above TeV scale but below the 
inflation scale, such as Pati—Salam unification. 
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